Solitons in coupled atomic-molecular Bose-Einstein condensates 
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We consider coupled atomic-molecular Bose-Einstein condensate system in a quasi-one- 
dimensional trap. In the vicinity of a Feshbach resonance the system can reveal parametric soliton- 
like behavior. We analyze bright soliton, soliton train and dark soliton solutions for the system in 
the trap and in the presence of the interactions between particles and find the range of the system 
parameters where the soliton states can be experimentally prepared and detected. 
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I. INTRODUCTION 



The experimental success in realizing a Bose-Einstein 
condensate (BEC) in dilute atomic gases [H, S| has 
led to the strong experimental and theoretical activity in 
this field. Dilute gases turn out to be an ideal system 
for realizing and high accuracy manipulation of various 
quantum many-body phenomena. The nonlinear nature 
of a BEC in the mean field description allows investiga- 
tions of e.g. vortex lattices [J, H, El , four- wave mixing 
0, Hj], Josephson-like oscillations U El and solitons in 
the system. The latter have been realized experimentally 
for a quasi one-dimensional (ID) atomic BEC with repul- 
sive effective interactions (so-called dark solitons) fill. 11211 
and with attractive interactions (bright solitons) [131 Il4j . 
There are also proposals to prepare experimentally soli- 
tons in a mixture of bosonic and fermionic atoms fl5j 
and in an atomic BEC with a time-dependent scattering 

length mm EE HI- 

In the present publication we consider another type of 
solitonic solutions. That is, parametric solitons, which in 
the context of nonlinear optics occur as a result of a cou- 
pling between fields in a nonlinear medium such that each 
field propagates as solitons [2l| . In the coupled atomic- 
molecular BEC, there is a nonlinear resonant transfer be- 
tween atoms and molecules, as well as terms proportional 
to the densities. The solitons of this type have been inves- 
tigated in nonlinear optics [22|, H3, 01 and in the problem 
of the self-localization of impurity atoms in a BEC [25| . 
In the context of the atomic-molecular BECs the para- 
metric solitons have been analyzed in Ref . [26|, [2tJ • The 
latter publications consider, in principle, solitons in free 
space in any dimension but they concentrate on the ex- 
istence and stability of the soliton like states in 3D free 
space. However, a profound analysis of the system in a 
quasi ID trap is missing and methods for experimental 
preparation and detection are not considered. Trapping 
an atomic BEC in a quasi ID potential allowed realiz- 
ing remarkable experiments where bright soliton wave 
packets travelled in ID space (still with transverse con- 
finement) without spreading [l3|, E]- In the case of the 
coupled atomic-molecular BECs, experimentalists have 
to face a problem of atomic losses, which becomes sig- 
nificant at a Feshbach resonance. We show that to over- 
come this obstacle one has to deal with moderate particle 



densities, and consider the range of the system parame- 
ters where preparation and detection of the solitons are 
sufficiently fast to be able to compete with the loss phe- 
nomenon. Moreover, the soliton states can be prepared at 
the magnetic field slightly off the resonance value, where 
the number of created molecules is still considerable but 
the rate of atomic losses much smaller than at the reso- 
nance. 

The paper is organized as follows. In Sec. II we intro- 
duce the model of the system and in Sec. Ill we analyze 
different solitonic states. Section IV is devoted to anal- 
ysis of methods for experimental preparation and detec- 
tion. 



II. THE MODEL 

In order to obtain a complete model of a Feshbach 
resonance in an atomic BEC, the intermediate bound 
states (molecules) have to be included explicitly in the 
Hamiltonian. At the resonance, the number of molecules 
becomes considerable and a second (molecular) BEC is 
formed. We take into account two body atom-atom, 
molecule-molecule and atom-molecule collisions, as well 
as the term responsible for the creation of molecules, i.e. 
transfer of pairs of atoms into molecules and vice versa 
HI: 
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where m is the atomic mass, A a , A m , A am are coupling 
constants for the respective interactions, and a deter- 
mines the strength of the resonance. The £ parameter 
is a difference between the bound state energy of two 
atoms and the energy of a free atom pair [281 ] . and it 
can be varied by means of a magnetic field. U a (f) and 
U m (r) stand for atomic and molecular trapping poten- 
tials, respectively. We would like to point out that in 
our model an internal structure of molecules and all pro- 
cesses involving the internal structure are neglected and 
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we consider only a single molecular state described by 
the operator ip m . 

In the present paper we consider the system param- 
eters mainly related to 87 Rb atoms. That is, the mass 
m, the coupling constant A a and the strength of the reso- 
nance a correspond to 87 Rb atoms and the Feshbach reso- 
nance that occurs at the magnetic field of B r = 685.43 G 
[29] ]. However, because the precise values of the cou- 
pling constants X m and X am are unknown these con- 
stants have to be chosen arbitrary. For simplicity we 



assume X a 



A m ~ X a — 47r?i a/m, where a = 5.7 nm 



is the value of the atomic scattering length far from a 
Feshbach resonance. In Sec. Ill we show that this as- 
sumption is not essential because in a broad range of 
X m and X am values there exist stable soliton solutions. 
The resonance strength parameter a can be expressed 
in terms of AB, the resonance width, and A/2, differ- 
ence between magnetic moments of a molecule and a free 



atom pair, a = y 4ttH aAjlAB/m. We have chosen for 
investigation the broad resonance which occurs at the 
magnetic field B T = 685.43 G where AB = 0.017 G and 
Afi = 1.4/x.b (mb is the Bohr magneton) [29| . We focus 
on the case where the system is prepared in a harmonic 
trap with so strong radial confinement that the radial fre- 
quency of the trap exceeds the chemical potential of the 
system. Then, only the ground states of the transverse 
degrees of freedom are relevant and the system becomes 
effectively ID. The chosen trap frequencies are u> mi ± — 
uj a j_ = 2ir x 1500 Hz and us m ^ x — us a ,x = 27r x 10 Hz 
[36j |. The effective ID coupling constants are obtained 
by assuming that atoms and molecules are in the ground 
states of the 2D harmonic trap of frequencies u> m .± and 
w a ,j_, and by integrating the energy density over trans- 
verse variables. In the harmonic oscillator units, 




(2) 



where Eq, xq and To are energy, length and time units, 
respectively, one obtains 



A Q — A,, 



0.505, 



and 



a 



41.0. 



(3) 



(4) 



The equations of motion for the atomic and molecular 
mean fields, corresponding to the Hamiltonian JT]), in the 
ID model, in the units ([2|), are 



dt 
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(5) 



where N is the total number of atoms in the system 
(28|. The ID detuning e is modified with respect to the 
corresponding £ in the 3D case, i.e. e = £ + {u} a ,± ~ 
2t^m,±)/^a,K = £ — 150. The wave-functions 4> a {x) and 
4> m {x) are normalized so that 



dx(|</. a (a;)| 2 + 2|0 m (x)| ; 



1. 



(G) 



The Eqs. ([5]) look like a pair of coupled Gross-Pitaevskii 
equations but with an extra term responsible for the 
transfer of atoms into molecules and molecules into 
atoms. The presence of this term allows for solutions 
that reveal parametric soliton like behavior. 

The mean field model (O can be used if the conden- 
sates are nearly perfect. That is, if the depletion effects 
are negligible which is usually true even for small par- 
ticle number in a system [30| . The model neglects also 
effects of particle losses. The latter can be described by 
introducing terms proportional to densities with imagi- 
nary coefficients (that can be estimated provided there 
is precise experimental analysis of the losses in the vicin- 
ity of the Feshbach resonance) in the square brackets of 

Eqs. © mLm. 



III. SOLITON LIKE SOLUTIONS 

A. Bright soliton solutions 

The time- independent version of Eqs. ([5]) reads 
1 d 2 1 



2dx 
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where /i is the chemical potential of the system, and we 
have assumed that the solutions are chosen to be real. 
The ground state of the system described by Eqs. ([7]) 
can be found numerically. However, in the absence of 
the trapping potentials and for A ? = Am_= X am = 0, 
there exists an analytical solution [22l. [23 . [24 1 if 



where 



18 \V3 



\a 2 NJ 



(8) 



(9) 



3 




'.5 -0.25 0.25 0.5 -3 -2 -1 1 2 3 

x x 



FIG. 1: Solid lines correspond to the ground state solutions 
of Eqs. © for TV = 100 (a) and 1000 (b). Dashed lines are 
the analytic solutions (|10p corresponding to the case without 
the interaction between particles and in the absence of the 
trapping potentials. Dotted lines are related to the gaussian 
ansatz (I13[) - for TV = 100 the states are almost the same as 
the exact ones. The atomic wavefunctions <f> a > 0, whereas 
the molecular ones <j>m < 0. The ground states of Eqs. Q cor- 
respond to the detuning chosen so that {<j> a \<j) a ) = (<t>m\<f>m), 
i.e. e = -900.53 (a) and e = -1505.25 (b). Note that the 
widths of the states are much smaller than the width of the 
ground state of the harmonic oscillator, which indicates that 
even for a very small particle number the nonlinearities in 
Eqs. determine the shapes of the states. 



That is, 



± 



cosh 2 (f ) 



cosh 2 (f ) 



(10) 



where 



(11) 



describe the parametric bright soliton of width I, which 
can propagate in time without changing its shape. The 



FIG. 2: Two bright soliton solution for TV = 100 (a) and 
TV = 1000 (b). The atomic wave-functions <j> a {x) are odd 
functions while the molecular wave-functions (f> m (x) axe the 
even ones. Note that populations of atomic and molecular 
condensates are unequal for the value of the detuning chosen, 
i.e. e = 0. Note that the widths of the states in panel (a) 
are much smaller than the width of the first excited state of 
the harmonic oscillator, which indicates that even for a very 
small particle number the nonlinearities in Eqs. ([7]) determine 
the shapes of the states. 



corresponding time-dependent solution is 



±- 



4>a(x,t) 
4>m(x,t) = - 



ae 



ivx p—i(v +^i)t 



cosh 2 (2^) 

ae 2ivx e -2i(v 2 + n)t 

cosh 2 : 



(12) 



where v is the propagation velocity. The analytical so- 
lution is known for a particular value of the detuning e 
but solutions revealing soliton like character exist also for 
other values. By changing e we can obtain solitons but 
with unequal population of atomic and molecular con- 
densates. 

In the presence of the (repulsive) interactions, the so- 
lutions for the ground states reveal wave-packets that, 
for increasing total number of atoms TV in the system, 
become wider and wider. Indeed, for A a > 0, A m > 
and X am > and with increasing TV, the repulsive in- 
teractions are stronger (see Fig. [TJ. Unlike the repulsive 
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nature of the interactions, the term responsible for the 
transfer between atoms and molecules acts as an attrac- 
tive potential. For N not too high, the transfer term 
dominates over the interaction terms and it alone de- 
termines the shape of the wave-packets. In Fig. [JJ we 
can see that the ground state solution for N = 100 is 
only slightly wider than the analytical solution (fTQ|) . In 
Fig. [1] we also present the results of variational analysis of 
Ref. (27| (including also nonzero coupling constants A m 
and A am ) where the solutions have been assumed to be 
given by gaussian functions, 



<j> a {x) = Ae~ 



4>b{x) = Be 



-bx 2 



(13) 



Requiring the same fraction of atoms and molecules, for 
N = 100 (N = 1000), we obtain A = 1.77, a = 139.54, 
B = 1.70, b = 117.43, s = -900 (I = 1.08, 5 = 19.58, 
B = 0.97, b = 12.47, e = -1429) which fits quite well to 
the exact solutions. 

Even in the presence of the interactions the obtained 
solutions do not loose their solitonic character, which 
can be verified by the time evolution of the states in 
ID free space (i.e. in the presence of the transverse 
confinement but with the axial trap turned off) — see 
Sec. IV for examples of the time evolution. The exis- 
tence of the solitons can be also verified employing the 
Gaussian ansatz (fT3"|) . Indeed, in the absence of the axial 
traps we have checked that, for N — 1000, A m S (0, 2A a ) 
and X am € (— 2A a , 2A a ), the solutions exist with a and b 
in the range (1,20000) [33j. It shows also that the choice 
of equal coupling constants is not essential in order to 
deal with soliton like solutions. 

All stationary states of the Gross-Pitaevskii equations 
(JT]) have been obtained numerically by means of the so- 
called imaginary time evolution which implies they are 
stable solutions of these non-linear coupled Schrddinger 
equations. 



B. Two bright soliton and dark soliton solutions 

In the absence of the trapping potential and for A a = 
A m = Aom = one can find an asymptotic solution of 
Eqs. (J7j that reveals two bright soliton structure, 
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FIG. 3: The soliton solution (solid lines), obtained for the 
detuning e — 0, compared to the solution within the Thomas- 
Fermi approximation (dashed lines). Panel (a) shows the 
atomic wavefunctions while panel (b) the molecular ones. The 
insets show the same as the main panels but with plots en- 
larged around the trap center. 



In the presence of the traps and the interactions we 
may look for similar solutions. Analysis of Eqs. (J7J shows 
that due to the fact that the trapping potentials are even 
functions we may expect that eigenstate functions are 
either even or odd. However, because of the presence of 
the </> 2 term in the second of Eqs. ([7]), an even function is 
the only possibility for <fi m . In Fig. [2] there are examples 
of the two bright soliton solutions for N — 100 and N — 
1000. 

With an increasing number of particles the shapes of 
the density profiles change. For large particle numbers 
the healing length of the system becomes much smaller 
than the size of the system and one obtains density 
shapes that resemble the dark soliton state known for 
the Gross-Pitaevskii equation in the case of repulsive in- 
terparticle interactions [HI, [13 — the results for N = 10 5 
are shown in Fig.[3J Here, however, only the atomic wave- 
function <p a (x) exhibits a phase flip. We can approximate 
the exact eigenstates by (see Fig. 
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are solutions for the ground state of the mean field equa- 
tions within the Thomas-Fermi approximation, and the 
chemical potential fi can be found from the normalization 
condition ©. 

Finally we would like to mention that in the case of the 
harmonic trap, having any solutions 4> a Q(x) and <j> m o(x) 
of the time-independent problem, Eqs. ([7]), correspond- 
ing to the chemical potential /i, one can easily obtain 
time evolution of the initial wave- functions 4> a o( x ~ <z) 
and 4> mQ (x - q), i.e. 



where 



and 



<f> a (x,t) 
<f>m(x, t) 



S(q) 



(b m0 (x-q)e^ t e 2 ^ x - s ^\ (19) 



to 



df [<f(t')-q 2 (t')], 



cfq 
dt 2 



q = 0. 



(20) 



(21) 



The proof can be done by direct substitution of (TIT))) into 
([5]). This indicates that, similarly as in the case of the 
Gross-Pitaevskii equation for a single condensate in a 
harmonic trap, time evolution of the translated station- 
ary solutions reveals harmonic oscillations of the center 
of mass of the particle cloud. 



IV. EXPERIMENTAL CREATION AND 
DETECTION 

To create experimentally a bright soliton state we pro- 
pose to start with a purely atomic condensate in a quasi- 
1D trap with the magnetic field significantly higher than 
the resonance value. Then the field should be slowly de- 
creased until it reaches the resonance. Actually, from the 
experimental point of view, it is better to end up slightly 
off the resonance because the rate of atomic losses is then 
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FIG. 4: Solid lines: density profiles of solitons excited after 
the magnetic field sweeping for N = 100 (a) and N = 1000 
(b). The sweeping lasted 95.5 ms (a) and 93.8 ms (b). Dashed 
lines show the exact ground states corresponding to the mag- 
netic field at the end of the sweeping. The detuning e — 
at the end of the evolution. Top curves, in each panel, are 
related to densities of atoms and bottom ones to densities of 
molecules. 



much smaller than the rate at the resonance. The final 
populations of atomic and molecular condensates are un- 
equal in such case but they still reveal solitonic charac- 
ter. In the following we analyze different velocities of the 
sweeping of the field in order to find the optimal shape 
of the state. 

We have done numerical calculations for N in the range 
between 100 and 10000. Assuming that the initial state 
is an almost purely atomic BEC, i.e. the ground state 
of the system far from the resonance, the detuning e has 
been varied linearly in time from 3331 to (Ae = 3331 
corresponds to AB = 0.017 G). In order to obtain a nice 
solitonic state, the evolution time has to be carefully cho- 
sen. We have found that the square overlap between a 
final state and a desired eigenstate oscillates as a function 
of the evolution time. Therefore one has to find a proper 
value of the evolution time, corresponding to a maximum 
value of the square overlap. From the experimental point 
of view, the shorter the evolution time the better because 
that allows one to reduce particle losses significant close 
to the resonance. In Fig. 0] there are results for N = 100 
and N = 1000 corresponding to the optimized evolution 
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FIG. 5: Soliton train obtained for N = 1000 after the sweep- 
ing of the magnetic field that lasted 25.5 ms. Solid line indi- 
cates the atomic density while the dashed line the molecular 
density. 



time and compared with the ground states at the reso- 
nance. One can see that despite quite short evolution 
times (95.5 ms and 93.9 ms for N = 100 and N = 1000, 
respectively), the final states reproduce the exact ground 
states quite well, i.e. the square overlaps are 0.94 and 
0.92. For N = 10000 the shortest evolution time that 
results in a reasonably high squared overlap (i.e. 0.83) 
is 90.7 ms. Note that the population of the molecular 
BEC is smaller than the population of the atomic BEC 
because the final value of the detuning corresponds to 
the magnetic field slightly off the resonance. 

Using the procedure described above it is possible not 
only to create a single bright soliton but also excite 
states that reveal multi-peak structure (so-called "soli- 
ton trains" [l4[ )• Actually, deviations from the optimal 
sweeping time of the magnetic field result in multi-peak 
structure of the atomic and molecular densities. An ex- 
ample of N = 1000, where the soliton train is particularly 
well reproduced, is shown in the Fig. [5] The calculations 
have been performed starting with the same initial state 
and the magnetic field value as in Fig. QJ but the sweeping 
time is now 25.5 ms. 

In Ref . [H [3 detection of time evolution of wave- 
packets without spreading in the presence of the trans- 
verse confinement but with the axial trap turned off (or 
even in the presence of an inverted axial potential) was an 
experimental signature of the bright soliton excitation in 
an atomic BEC with attractive atom-atom interaction. 
In the present case we can use the same test to check 
if the states obtained reveal stable soliton like charac- 
ter. It turns out that indeed all the states presented do 
not spread when the axial trap is turned off during time 
evolution lasting even 100 ms. We do not attach the cor- 
responding figures because they would show practically 
the same density profiles as the densities of the initial 
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FIG. 6: Time evolution of a soliton state for N = 1000 in 
an inhomogeneous magnetic field (B gra d = 1 G/cm), for the 
detuning (I25[) e gra d = 0. Panel (a) shows the atomic wave- 
packet, panel (b) the molecular one. The time t = 1 in the 
units ([2]) corresponds to 15.9 ms. 



states. 

In Ref. [34{ the results of experimental creation of 
molecules from an atomic BEC prepared not in a quasi 
ID potential but in a 3D trap are presented. After 
producing the molecules, the experimentalists analyze 
the evolution of the system with the traps turned off 
but in the presence of an inhomogeneous magnetic field. 
They observe separation of the molecular cloud from the 
atomic one, which is due to the difference between atomic 
and molecular magnetic moments. We will analyze the 
evolution of the atomic and molecular clouds in an in- 
homogeneous magnetic field but in the presence of the 
transverse confinement. One may expect that when the 
solitonic states are created there is strong attractive cou- 
pling between atoms and molecules that will compete 
with the opposite effect resulting from the difference of 
the magnetic moments. 

In the magnetic field gradient, eqs. ([5]) are modified. 
Instead of the detuning e there are terms proportional to 
the magnetic moments: p> a {B gra( iXoX + B)/Eq - in the 
first equation, and fx m {B gra( iXQX + B) / E$ - in the second 
one. -Bgrad is a value of the magnetic field gradient and 
xq and Eq are length and energy units ©. By applying 
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FIG. 7: Time evolution of a soliton state for AT = 100 in 
an inhomogeneous magnetic field (B gra d = 1 G/cm), for the 
detuning (|25[) £ gra d = 2216. Panel (a) shows the atomic wave- 
packet, panel (b) the molecular one. The time t = 1 in the 
units (j2} corresponds to 15.9 ms. 



FIG. 8: Time evolution of a soliton state for N = 1000 (b) 
in an inhomogeneous magnetic field (B gra d = 1 G/cm), for 
the detuning (|25[) £ gra d = 2216. Panel (a) shows the atomic 
wave-packet, panel (b) the molecular one. The time t = 1 in 
the units ([2]) corresponds to 15.9 ms. 



suitable unitary transformations, however, these terms 
can be eliminated from the first equation. First, the co- 
ordinate transformation x — > x + f3 a is applied, then the 
phases are adjusted according to: 

^ a -> Vae* ( ^"^ )4 (22) 

Vn, -» ^-^ ) ', (23) 

where /3 a = x<jjl a B g!:a d/ Eq. Finally, there are only two 
additional terms in ([5]) (both present in the second of 
eqs. ©): 

(3x = ^ Bgrad x, (24) 

and a detuning: 

£grad = £ - f -|r-Bgrad^ j3 a A/l, (25) 

modified with respect to the case without gradient (i.e. 
with only e). 

The behavior of the system in the presence of an in- 
homogeneous magnetic field depends on the value of the 



detuning e gra d. For e gra d = the clouds do not separate 
but there is an efficient transfer of atoms into molecules. 
In Fig. M we see that for N = 1000, S grad = 1 G/cm 
and Egrad = 0, the number of molecules is increased af- 
ter 15.9 ms of time evolution. For e gra d = 2216 (cor- 
responding to e = 0) and for small particle number 
(e.g. N — 100) there is an opposite transfer, that is 
of molecules into atoms (see Fig. [7]), and when molecules 
are absent the atomic wavepacket starts spreading. How- 
ever, if N is greater (N = 1000 or more) the wavepackets 
begin to split into several separated peaks (see Fig. [5]) but 
there is still no separation between atomic and molecu- 
lar clouds. Similar splitting of solitonic wavepackets has 
been analyzed in the case of a single component BEC in 
the presence of a gravitational field [35] . For greater field 
gradient than the value chosen in Figs.[6][Sor slightly dif- 
ferent values of the coupling constants we observe quali- 
tatively similar behaviour. 

The longer creation and detection of the solitons lasts, 
the more serious becomes the problem of atomic losses in 
the vicinity of the Feshbach resonance. Measurements of 
the losses for different Feshbach resonances in 87 Rb were 
done in a 3D trap (frequencies: 2-7r x 50 Hz, 2-7T x 120 Hz, 
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2ir x 170 Hz) with an initial total number of atoms 
N = 4 x 10 6 (that corresponds to the density at the cen- 
ter of the atomic condensate of order of 10 15 cm -3 ) [29l |. 
The observed loss rate, defined as a fraction of atoms 
lost during 50 ms hold time in the trap, was 78% for 
the resonance at B r — 685.43 G. In the present publica- 
tion we have chosen very high frequency of the transverse 
traps (i.e. 2ir x 1500 Hz) in order to ensure that the ID 
approximation is valid for all values of particle numbers 
considered. It results in the peak density for solitonic 
states comparable to the experimental density in [2^ |. 
However, the density can be significantly reduced if one 
chooses weaker transverse confinement and a moderate 
particle number. For example, for N = 1000 and for the 
transverse traps of 2tt x 200 Hz, the peak soliton density 
is 5 x 10 14 cm" 3 and the evolution time needed to obtain 
the final state with 0.91 squared overlap with the desired 
eigenstate is 74.8 ms. The losses will be smaller if the 
magnetic field at the end of the sweeping corresponds to 
a value slightly above the resonance, where the solitons 
still form but the loss rate is smaller. 

Actually, solitons can be prepared and detected in the 
presence of the losses provided that the escaping atoms 
do not significantly excite the system of the remaining 
particles. Indeed, an experimental signature of the soli- 
ton excitation could be the time evolution of the system 
with the axial trap turned off, where the wavepackets 
are expected not to increase their widths. Importantly 
when TV decreases, the influence of the atom-molecule 
transfer term (which is responsible for the existence of 
the solitonic solutions) on the character of the solutions 
grows, as compared to the influence of the interaction 
terms. This is because the former depends on the par- 
ticle number as y/N while the latter as N, see Eq. ([5]). 
Consequently the normalized particle density should not 
increase its width even in the presence of the losses if a 
soliton state is prepared in an experiment. 

V. SUMMARY 

We consider solitonic behaviour of coupled atomic- 
molecular Bosc-Einstein condensates trapped in a quasi- 
1D potential. In the vicinity of a Feshbach resonance 
there are eigenstates of the system that reveal soliton- 
like character. The profiles of the solitons depend on the 



total number TV of the particles in the system and on the 
values of the scattering lengths characterizing particle in- 
teractions. For moderate particle numbers and positive 
scattering lengths one can find two bright soliton solu- 
tions which for large N turn into a state that resembles 
a dark soliton profile known for atomic condensates with 
the repulsive atom-atom interaction [111, [l2j • 

We analyze methods for experimental preparation and 
detection of the solitons in the atomic-molecular BECs. 
They can be obtained experimentally starting from a 
purely atomic BEC and slowly decreasing the magnetic 
field down to a value slightly above the resonance. Evo- 
lution times necessary to create solitons are quite short 
which is very promising experimentally because it allows 
reducing atomic losses important in the vicinity of the 
resonance. To detect the solitons one can perform (sim- 
ilarly as in Refs. [H, [IH) the evolution of the system in 
the absence of the axial trap (but still in the presence of 
the transverse confinement) where the wave-packets are 
expected to propagate without spreading. 

We have also considered creation of the bright soliton 
trains in a coupled atomic-molecular condensate. Soli- 
ton trains in Bose-Fermi mixtures obtained by temporal 
variations of the interactions between the two compo- 
nents has been recently proposed Our system is 
different but its qualitative behaviour is to some extent 
similar. In both cases, there is a term in the Hamil- 
tonian that leads to an effective attractive interactions. 
This term can dominate over repulsive interactions and 
stabilize bright solitons. 
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